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ABSTRACT 


In  this  paper  we  study  a  model  for  the  population  transition  probabilities 
for  a  branching  process  composed  of  particles  diffusing  in  a  finite  interval.  The 
model  is  in  general  non -Markovian  since  we  assume  the  branching  transformation 
probabilities  (for  a  particle  depend  on  its  age  and  position.  The  process  is 
described  by  the  random  number  NJx)  of  particles  in  the  interval  I  at  time  t 
that  are  generated  by  a  single  particle  initially  at  the  point  x  in  I.  By 
considering  N^x)  as  a  regenerative  process  with  respect  to  the  random  age  and 
position  of  the  initial  particle  when  it  is  transformed,  we  develop  a  functional 
equation  for  the  generating  function  for  NJx)  .  This  functional  equation  is  the 
basis  for  our  study  of  the  population  probabilities  P(Nt(x)  =  n)  ,  n  =  0,  1,2,..., 
as  function  of  x  in  I  and  t  in  [0,  00 )  . 

Our  principal  result  concerns  the  extension  to  our  model  of  the  fundamental 
result  on  the  probability  for  ultimate  extinction  in  a  Galton -Watson  branching 
process,  [3,  Chapter  XII]  i  Letting  f Q(x)  be  the  probability  for  ultimate  extinction 
in  a  process  generated  by  a  particle  initially  at  x,  we  give  a  necessary  and 
sufficient  condition  for  fQ{ x)  =  1 .  We  also  give  a  characterization  of  fQ(x)  , 
x  in  I,  as  the  minimal  positive  solution  of  a  functional  equation  which  is 
asymptotically,  related  to  the  basic  equation. 


EXTINCTION  PROBABILITIES  FOR  AGE  AND  POSITION -DEPENDENT  BRANCHING  PROCESSES 


Howard  E.  Conner 

Introduction  . 

The  process  is  composed  of  particles  diffusing  in  a  bounded  open  interval 
1^,  and  it  is  generated  by  one  particle  initially  at  x  in  I  .  The  diffusion  is 
represented  by  the  conditional  probability  density  function  (  p.  d,  f)  ■  p(x,  y;  s)  ,  • 
where  p(x,  y;  s)dy  is  the  probability  for  a  particle  at  x  at  time  t  to  be  at  y 
at  time  t  +  s,  assuming  it  has  not  been  previously  to  the  boundary  D  of  IQ  . 
Therefore,  we  assume 

1.1  (a)  p(x,  y ;  s)  =  0  for  x  in  F,  y  in  I  and  0  <  s  <  <*>  and 

(b)  p(x,  y;  s)  is  non-negative  and  continuous  for  x,  y  in  I  and 
0  <  s  <  °°  . 

The  boundary  T  is  an  absorbing  barrier  for  which  a(x,  t)  is  the  probability 

distribution  function  for  a  particle  initially  at  x  to  be  absorbed  at  r  within  the 

time  t  and  a(x)  =  lim  a(x,  t)  is  the  probability  for  the  ultimate  absorption  at  T 
t-*oo  . 

for  a  particle  at  x.  Consequently,  we  assume 

Sponsored  by  the  Mathematics  Research  Center,  U.  S.  Army,  Madison,  Wisconsin,- 
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1.2  (a)  a(x,  t)  is  non-negative  and  continuous  for  x  in  I  and 

0  <  t  <  °°  with  a(x,  t)  =  1  for  x  in  r  and  0  <  t  <  00 , 

( b)  for  each  x  in  I,  a(x,  t)  is  nondecreasing  as  t  increases  and 

(c)  a(x)  =  lim  a(x,  t)  is  continuous  for  x  in  I  with 

t-*-00 

a(x)  =  1  for  x  in  r  and.  not  identically  1  for  x  in  I  .  . 

Each  particle  has  an  independent  random  life  span  with  the  p.  d.  f.  g(  s)  where 
g(  s)  ds  is  the  probability  for  a  particle  formed  at  time  T  to  end  its  life  at  time 
T  +  s  ;  and  so,  we  assume 

1.3  g(  s)  is  non-negative  and  continuous  for  0  <  s  <  °°  . 

A  particle  whose  life  has  ended  at  the  age  T  and  at  the  point  x  in  IQ  is  replaced 

by  a  random  number  k  of  particles  with  probability  b,  (x,  T)  and  the  expected 

co  k 

number  of  particles  with  which  it  is  replaced  is  Yj  kb  (x,  t)  .  Therefore.,  we 

k=l  * 

assume 

(a)  b  (x,  s)  is  a  non-negative  continuous  function  for '  x  in  I  and 

K  U 

.0  <  s  <  , 


1.4 
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00 

(b)  Ij  b  (x,  s)  =  1  for  x  in  I  ,  0  <  s  <  00 . 
k=  0  K  0 


00 

(c)  [i(x  ,  s)  =  l  kb  (x ,  s)  is  positive,  bounded  and  continuous 
k=  1 


for  x  in  I  and  0  <  s  <  00 . 


At  any  given  time  t  a  particle  initially  at  x  in  I 


1.5  (a)  has  been  absorbed  at'  T  with  probability  a(x,  t)  or 


(b)  has  ended  its  life  in  I  with  probability  f  f  p(  x,  y;  s)  g(  s)  dsdy  or 

'  I  0 


( c)  is  diffusing  in  IQ  with  probability  b(x,  t)  . 


Ultimately  a  particle  is  either  absorbed  at  T  or  transformed  into  new  particles  at 
some  interior  point  x.  Consequently,  we  assume 


1.6  b(x,  t)  is  non-negative  and  continuous  with  lim  b(x,  t)  =  0  as  t-+  00 
for  x  in  I  and  0  <  t  < 00 . 


Since  the  events  described  in  1.5  are  mutually  exclusive  and  exhaustive,  we  assume 


t  ... 

(a)  ■  1  =  a'(x,  t)  +b(x,t)  +  f  f  p(x,  y;.  s)  g(  s)  dsdy  and 

10  . 


1.7 
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oo  . 

(b)  1  =  a(x)  +  ff  p(x,  y;  s)g(s)dsdy 

I  0 

for  x  in  I  and  0  <  t  <°°,  where  the  integral  in  (b)  is  assumed  to  converge  ■ 
uniformly  for  x  in  I . 

We  define  the  generating  function  for  the  transition  probabilities 

fk(x,  t)  =  P(Nt(x)  k)  ,  k  =  0,  1,  2,  ...  ,  by 
00 

1.8  f(x,  t;  z)  =  Yj  Mx>  t)2^  >  H'x,  t;  1)  =  1, 

k“  0 

defined  for  x  in  I,  0  <  t  and  lz|  <1.  if  the  initial  particle  is  replaced  at 

the  age  s  by  k  particles  at  the  point  y  in  IQ,  the  generating  function  for  the 

k' 

population  size  at  t  >  s  is  then  f  (y,  t  -  s  ;  z)  .  Since  the  probability  for  the 
initial  particle  to  be  replaced  at  an  age  s  with  k  particles  at  y  is  given  by 
b  (y,  s)p(x,  y;  s)g(s)dsdy,  this  and  the  previous  discussion  suggests  f(x,  t;  z) 
satisfies  the  functional  equation, 

t 

1.  9  f(x,  t ;  z)  =  a{x,  t)  +  zb(x,  t)  +  //  h[y,  s;f(y,  t-s;  z)]p(x,  y;  s)g(s)dsdy 

10 

for  x  in  I,  0  <  t  <  00  and  [  z [  <  1,  where 

oo 

1.  10  h(x,  t;  z)  =  Yi  t)  zk  I  h(x»  t;  1 )  =  1 

k=0  K  . 

Having  formally  derived  1.9,  it  is  necessary  to  establish  the  existence  of  a 
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. unique  solution  f(x,  t;  z)  bounded  by.  1  which  is  a  generating  function  in  z 
for  each  x  and  t. 

Before  doing  this,  we  introduce- 

;°o  ’  ■ 

1.11  K(x,  y )=  /  H-(y,  s)p(x,y;  s)g(  s)ds, 

the  expected  number  of  particles'  replacing,  a  particle  initially  at  x  which  ends; 
its  life  at  y.  We  infer  from  the  conditions  of  nonnegativity,  continuity  and  uniform 
integrability  that  K(x,  y)  is  nonnegative  and  continuous' for  x,  y  in  I.  Therefore. 
K(x,  y)  defines  a  linear  integral  operator  K 

1.  IE  Kf(x)  =  f  K(x,  y)f(y)dy  ,  x  in  I 

I 

(a)  transforming  the  class  of  functions  continuous  on  I  into  the 
subclass  of  functions  vanishing  on  T-  and 

( b)  transforming  the  convex  cone  of  nonnegativ.e  functions'  on 
I  into  itself. 

Consequently,  the  operator  K  has  special  spectral  properties  .similar  to.  those  for 
a  matrix  operator  with  nonnegative  elements.  Of  these  we'use  .the  existence  of 'a 
positive  characteristic  number  of  minimum  modulus  to  form  .a  necessary 'and' sufficient 
condition  for.  ultimate  extinction  for  Nt( x)  with- probability.  1 . 
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A  number  of  mathematical  models  for  branching  processes'  are  covered 
in  the  papers  of  T.  E.  Harris-,  [4] .  -Limiting  theorems  for  specific  models  with' 
transformation  probabilities  depending  on  age' or  position  are  also  given' in  these 
papers'and  in  papers  by  R.  Bellman  and'T.  E.  Harris,-  [1],  N.  Levinson,  [6]  and 
the  author;  [  2] .  The  problem  of  ultimate  extinction  for  a  time -dependent  branching 
process  is  given  in  papers  by  B.  A. •  Sevast'janov,  [  7,  8].  '  ■ 

§2.-  Existence  of  generating  function.  ' 

The  existence  of  a  solution  to  1.  9  is  given  by  the 

Theorem  2.  1.  Suppose  a(x,  t),  b(x,  t),  p(x,  y;  t) ,  g(t),  h(x/t;  z)  and -  |i(x,  t) 
satisfy  the  conditions  given  in  1.  1-1.  4,  1.  6,  1'.  7  and  1.  10-.  With  these  conditions, 
there  exists  a  unique  solution  f(x,  t;  z)  to  1.  9  defined  for  x  in  I/  O  <  t  <  00  ■  ' 
and  I  zl  _<  1  and  satisfying: 

"(a)  f(x,t;.z)  is  continuous  on  its  domain, 

(b)  lf(x,  t;  z)  I  <1  with  f(x,  t;  1)  =  1  on  its  domain  and 

(c) ,  for  each  x  in  I  and  t  >  0,  f(x,  t;  z)  is  analytic,  on 

I  z|  <  1  and  therefore  has  a  representation 

f ( x,  t;  z)  =  Yi  f^x>  ^  zk 
k=0 

where  each  f  (x,  t)  is  continuous  and  nonnegative  for  x  in  I . 

■  JC 

'  and  0  <  t .  ■ 
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Proof.  We  first  list  some  properties  of  h(x,  t;  z)  which  follow  directly  from  ■ 
the  condition  1.4  and  1. 10  and  an  application  of  the  mean  value  theorem: 

2.1  (i)  0  <  h(x,  t ;  t^)  .<  h(x,  t ;  t^)  , 

(ii)  0  <  I  h( x ,  t ;  z)  |  '<  h(x,'t;  i  z I),  . 

(  iii)  0  <■  I h(x,  t ;  z)  1  <  .|x(x,  t).  .  and  .  ‘  „ 

(  iv)  |h(x,  t ;  z)  -  h(x,  t;  w)  |  •  <  |i(x,  t  )  I  z  -  w| 
for  x  in  I,  0  <  t,  I  z|  and  I  w[  <  1  and  0  <  -r  <  t'  <  1  . 

X  .•  La 

We  define  the  sequence  {f^Cx,  if  z)  }  with  fQ(x,  t;  z)  =  0  by 

t 

2.2  f  ,  (x,t;z>=  a(x,t)  +  zb(x,tj  +/  /  h[  y,  s;  f  (y,  t  -  s,  z)  ]  p(x,  y;.  s)  g(  s)  dsdy 
nrl  ,10 

for  x  in  I,  0  <  t  .and  |z|.<  1. 

An  induction  argument  using  properties  2.1  and  1.7  implies  for 
each  n ,  n  =  0 ,  1,2,..., 

2.3  .  |  fn(  x ,  t ;  z)  1  <  1  •  and  fn(x>  t;  1)  =  1 
for  x  in  I,  0.<  t  and  |z|'  <  i  . 
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A  second  induction  argument  using  the  continuity  properties  of  the  given, 
functions- and  the  analyticity  of  h(  x,  t  ;  z)  in  z  shows  each  f  (x,  t;  z)  is 
continuous  for  x  in  I,  t>  0  and  I  z[  <  1  .  and  for  each  x  in  I  and'' t  >  0 
it  is  analytic  for  |z|  <1.  Another  induction  argument  using  the  nonnegativity 

a  k 

of  the  given  functions  shows  the  derivatives  )  f^(x,  t;  0).  are  continuous 
and  nonnegative  for  k  =  0,  1 ,  2,  ...  .  Therefore  we  assert  each  f^(x,  t;  z) 
has  the' representation  for  x  in  I,  0  <  t  and  I  z|  <  1  given  by 

00 

V'  k 

2.4  f  (x,t;z)  =  /j  f  ,  (x,  t)z  with 

n  k  -  o  n*  K 

f  ,  (x,  t)  continuous  and  nonnegative,  k=0,l,2,...  .. 
n,  k 

We  use  2. 1  (  iv)  and  2.2  to  make  the  following  estimate,  for  general  n 


L 

^n+/X,t;Z)“fn(X,t;Z)  —  //^n^»t"s;z)''fn_i^>t-s>z)tlJ-(  y>s)P(x>y;s)g(s)dsdy 


on  x  in  I,  0  <  t  and  |z|  <  1.  The  substitution  of  2.3  with  n  =  1  into  the 
integrand  of  2. 5  with  n  =  1  gives 


t 

If  (x,  t;  z)  -  f  (x,t;z)l  </  / ^(y,s)p(x,  y;  s)g(s)dsdy. 

0  1 


Setting  A(T)  =  max  {  / h-(  y,  s)  p(x,  y,  s)  g(  s)  dy  I  x  in  I,  0  <  s  <  T  }  for  0  <  T, 

I 


we  have  If  (x,  t;  z)  f  (x,  t;  e)  1  <  t A  (  x)  ;  and  so,  by  induction,  we  have  for 

n  -  0,  1 )  2>^  •••  ^ 

2.6  I  fn  +  i  ( x  *  t »  z)  -fn(x,t;;z)|  <  tnAn(T)/n!' 

on  x  in  1 ,  0  <  t  <  T  and  f  z  |  <  1 , 

We  infer  from  2.  6  that  the  fn(x»  t;  z)  converge  as  n-r0®  to  a  limit 
f(x,  t;  z)  and  that  this  convergence  is  uniform  for  any  T  >  0  on  x  in  '1,0  <  t  <  T 
and  |  zf  <  1  .  Therefore  using  the  Moore -Weiers trass  theorems,  we  know 

(a)  f(x,  t;  z)  is  continuous  for  x  in.  I,  t  <  0  and  |'z[  <  1, 

(b)  f(x,  t;  z)  is  analytic  in  |z|  <1  for  each  x  in  I,  0  <  t  and. 

by  letting  n^10  in  2,2, 

(c)  f(x,  t ;  z)  is  a  solution  to  1.  9  for  x  in  I,  0  <  t  and  I'zi  <  1. 

The  properties  [f(x,  t;  z)  I  <  1  and  f(x,  t;  1)  s  1  are  consequences  of  the 

convergence  and  2.3.  We  use  the  Cauchy  formula  to  write  for  each  x  in  I  and 
0  <  t, 

2tt 

2.7  fn,k(X}  t)  “  ^  fn(x»  t?  el0)e"lk0d0* 

Since  fn(x>  t ;  z)  converges  uniformly  on  x  in  I,  0  <  t  <  T  and  |z|  <  1  for  each 
T  >  0,  we  know  the  right  side  of  2.7  converges  uniformly  on  the  same  sets  to 


vx’  l> a  k? 


'4\ 


a^j  f(x,  t;  0)  .  Therefore,  for  any  T' >  0  each  f^(x,  t)  is' the 

uniform  limit  on  x  in  I  ,  O' <  t  _<  T  and  I  z|  <  1  of  the  continuous  nonnegative 

functions  f  ,  (x,  t) .  given  in  2. 4  ;  and  so,  we  know  each  f  ix,  t)  is  ' 
n,  k  ■  .  k 

continuous  and  nonnegative  for  x'  in  I  and  0  <  t. 


To  complete  the  proof  we  must  show  the  solution  f(x,  t;  z)  is  unique 
in  the  class  of; functions  continuous  and  bounded  in  magnitude  by  1  on  x  in  I, 
0  <  t  and  i  z |  <  1  .  Suppose  g(x,  t ;  z)  is  a  continuous  function  with 
lg(x,  t;  z)|  <  1  on  x  in  I,  0  . <  t  and  I  z l  _<  1  satisfying  1.9  and  consider 
r( x,  t ;  z)  =  |  f(  x,  t ;  z)  -  g(x,  t ;  z)  |  .  Since  I  g( x,  t ;  z)  |  <  1  we  can  use 
1.9  and  2.  l(iv)  to  make  the  estimate 


t 

r(x,  t,  z )  <  /  /  r(  y,  t-s,  z)|o.(y,  s)p(x,  y;  s)g(s)dyds. 
0  I 


Letting  R(T,  X,  z)  =  max  {r(  x,  t ;  z)  e"  Xt  I  x  in  I;  0<t<T}  for 
T  and  ^  >  0  and  |z|  <  1,  we  have  R(T,  X,  z)  <  °o  and 


00 

R(  T,X,  z)  <  R(T,  X,  z)  /  /e  XS(j.(y,  s)p(x,  y;  s)g(s)dyds. 

0  1 

We  know  the  integral  exists,  is  finite,  1.4  and  1.7,  and  is  a  decreasing  function 
of  X;  so,  we  choose  XQ.  so  that  the  value  of  the  integral  is  less  than  j  .  This 
shows  0  _<  R( T ,  Xq ,  z)  _<  0  for  each'  T  >  0  and  |z|  <  1  and  consequently 
g(x,  t;  z)  =  f(x,  t;  z)  .  Therefore  f(x,  t;  z)  is  a  unique  continuous  solution 
to  1.9  satisfying  |  f(  x,  t ;  z)  |  <  1 ,  completing  the  proof. 
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Similar  functional  equations  can  be  derived  and  solved  for  the  multivariate 
generating  functions 

k  k. 

F(x,  tx,  ,  tn;  Zj,  .i.  ,  zn)=  Z.  P(Nt  (x)  =  k^;...  ;Nt  (x)  =  k^)zl  ...  z^ 

kj,  ...  *  k.  1  n 

satisfying  the  consistancy  condition  .that  for  every  set  of  positive  integers 

(  i .... ,  i  )  ,  1  <  m  <  n, 

1  m  —  — 

T(x,  tj,  ...  ,  tn;  zL,  ...,  Zn  =  F(x,  tv  ...  ,  tm;z  ...,z,  ) 

1  m 

where  z.  =  1  if  i  is  not  in  (i, ,  ...  ,  i  )  . 
i  1  '  m 

§  3.  Probability  for  eventual  extinction. 

Having  established  the  existence  of  a  unique  continuous  solution  f(x ,  t;  z) 
to  1.9  on  x  in  I,  0  <  t  and  |z|  <1  satisfying 

00  00 

3.1  f(x,.t;  z)  =  Z  fJx>  t)  z k  ;  Z  f  (x,  t)  =  1 

k=0  •  k=0 

where' each  f^(  x.  t)  is  nonnegative  and  continuous,  we  formally  identify  the 

f  (x,  t)  as  the  population  transition  functions  for  the  population  size  N  .(x) 
x  t 

for  a  system  generated  by  a  particle  initially  at  x .  In  particular  setting  z  =  0 


in  1 . 9  ,  we  have 
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t. 

3.2  f  (x,  t)  =  a(x,  t)  .+  ffh[y,  s;f  ( y,  t  -  s)]p(x,  y;  s)g(  s)dyds  . 

1  0 

where  fQ(x,  t)  formally  gives  the  Prob.,  {N^x)  =  o|n  (x)  -  1}.  Since  fQ(x,  t) 
is  nonnegative,  bounded  above  by  1  and  continuous  on  x  in  I  and  0  <  t, 

Theorem  2.1,  we  know  from  the  uniqueness  argument  following  2.8  that  fQ(x,  t) 
is  a  unique  continuous  solution  to  3.2  satisfying  |fQ(x,  t)  I  <  1  .  We  will  use 
this  uniqueness  property  to  show  fQ(x,  t)  is  an  increasing  function  of  t  for 
each  x  in  I,  a  property  which  is  suggested  by  the  probability  context. 

We  define  a  sequence  of  functions  for  x  in  I  and  0  <  t  by  g  (x,  t)  =  0 
and  g^+^(x,  t)  =  Ttj^x,  t)  where  T  is  a  non-linear  operator  defined  by 

t 

3.3  Tg(x,t)  =a(x,  t)  +  /{  /h[  y,  s;  g(  y,  t  -  s)  ]  p(x,  y;  s)  g(s)  ds}dy, 

I  0 

for  a  continuous  function  g(x,  t)  on  x  in  I  and  0  <  t  . 

An  induction  argument  shows  each  9n(x>  t)  is  a  nonnegative  continuous  function 
bounded  above  by  1  on  x  in  I  and  0  <  t,  and  it  is  a  monotone  increasing  function 
of  t  for  each  x  in  I.  The  last  statement  follows  directly  from  the  nonnegativity 
of  the  given  functions  and  the  nondecreasing  behavior  of  a(x,  t)  as  a  function 
of  t,  1.2.  With  the  same  estimates  as  used  in  Theorem  2.1  we  have 
[g1(x,  t)  -  gQ(x,  t)  [  <  1  and 
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I  g  (x,  t)  -  g  (x,  t)  [  <  /  /  [g  (y,  t-s)  -g  ( y,  t-s) I  \l(  y,  s)  p(x,  y;  s)g(  s)  dyds  . 
rvi  -i  n  X  0  ^  ^ 

Therefore  we  have  by  induction  on  n  ;  for  each  T  >  0  ,|  g  (x,  t)  -g  (x,  t)|<  tnAn(T)  /n 

for  x  in  I  and  0<  t<T  where  A(T)  is  defined  in  2.6.  Using  the  uniqueness 
of  fQ(x,  t)  as  a  continuous  bounded  solution  to  3.2,  this  shows  fQ(x,  t)  is  for 
each  T  >  0  the  uniform  limit  as  n  -*•  00  for  x  in  I  and  0  <  t  <  T  of  the  sequence 
gn(x,  t)  .  Consequently,  fo(x,t)  is  a  monotone  increasing  function  of  t  for  each 
x  in  I .  Summarizing  we  have 

3.4  ^  i's  a  nonne9ative  continuous  unique  solution  to  3.2  bounded 

in  magnitude  by  1  and  for  each  x  in  I,  it  is  a  monotone  increasing 
function  of  t. 

In  particular  3.4  implies  ^lim^  f Q(x,  t)  =  fQ(x)  exists  for  each  x  in  I. 

The  probability  context  suggests  fQ(x)  is  the  probability  for  eventual  extinction  in 

a  process  generated  by  a  particle  initially  at  x.  If  we  can  show  f  (x,  t)  is  continuous 

in  x  uniformly  in  t  then  we  know  f  (x)  is  continuous  in  x  and  therefore  by 

0 

the  Moore-Osgood  theorem  on  iterated  limits 

3.  5  lim  f  (x,  t)  =  f  (x)  <  1 

t  00  U  '  U 

uniformly  for  x  in  I .  To  show  f  (x,  t)  is  continuous  in  x  uniformly  in  t ,  we 


make  the  estimate 
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.•  "  t  •  '  ;  ■ 

If^t)  -f0(x2>t^f  latent)  -a(x2,t)l  +  //h[y,s;f0(y,t-s)]|p(x1,y;s)-p(x2,y;s)|g(s)dyds 

00 

<  |a(x  ,t)  -  a(x  ,t)  |  +  /  /  I  p(x  ,y;  s)  -  p(x  ,y;s)  |  g(  s)  dyds  . 

1  0  I  1  ' 

The  uniform  integrability  condition  1.  7  gives  the  existence  of  a  positive 
■y  { e ) ,  0  <  <=  ,  such  that  the  2nd  term  is  <  if  I  x  -  x  I  <y  )  .  Since  we 

u  i-  L,  Ct 

have  assumed  lim  a(x,t)  =  a(x)  uniformly  in  x,  1.3,  there  exists  a  positive 

t-»00 

6(e),  0  <  e  ,  such  that  the  first  term  is  if  Ix^-  xj  <  6(^-)  .  This  completes 

the  argument  and  establishes  3.5. 

Letting  t  -*■  in  3.2  and  using  1.7  and  3.5  shows  fQ(x)  is  a  continuous 
solution  to 

00 

3.  6  fn(x)  =  a(x)  +  /{  /  h[y,  s;f  (y)]p(x,  y;  s)g(  s)ds}  dy 

10 

with  0  <fQ(x)  <1  for  x  in  I.  Since  the  function  l(x),  l(x)  =  1  for  x  in  I, 
is  a  solution  to  3.  6  by  1.7,  we  have  the  problem  of  finding  a  necessary  and 
sufficient  condition  for  3.6  to  have  a  nonnegative  continuous  solution  different 
from  l(x)  . 

We  define  the  nonlinear  Urysohn  operator  U  on  the  class  of  nonnegative 
continuous  functions  bounded  in' magnitude  by  1  on  I, 
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CO 

3.7  ■  Uf(x)=a(x)  +  /  {/  h[y,  s;  f(  y)  ]p(x) ,  y;  s)  g(  s)  ds}dy . 

10 

The  proof  of  the  next  Theorem  is  based  on  the  relationship  between  U  and  the 
nonlinear  operator  T  defined  in  3.  3  and  the  spectral  properties  of  the  derivative 
of  .  U  at  1(  x)  =  1 ,  which  is  defined  by 

00 

Kf(x)  =  /{  /  h(y,  s;  1)  p(x,  y;  s)  g(  s)  ds  }f(  y)  dy  =  /k(x,  y)  f(  y)  dy , 

I  0  I 

where  K(x,  y)  is  the  expected  number  of  particles  replacing  at  y  a  particle  initially 
at  x.  The  transformation  properties  of  K  are  given  in  1.12,  and, as  stated 
there,  K  has  some  special  spectral  properties. 

The  spectral  properties  we  want  can  be  most  easily  obtained  by  approximating 
K(x,  y)  and  the  continuous  functions  on  which  it  operates  by  step  functions, 
forming  a  matrix  operator  with  nonnegative  elements  to  which  the  Perron-Frobenius1 
theory  can  be  applied.  A  more  general  theory  is  developed  in  the  monograph  by 
M.  G.  Krein  and  M.  A.  Rutman,  [5].  We  now  list  for  convenience  those  properties 
for  which  we  have  use.  Specifically,  assuming  K(x,  y)  is  nonnegative,  not 
identically  zero  and  continuous  for  x,  y  in  I,  the  following  statements  are  valid. 

3.8  The  kernel  K( x,  y) has  a  positive  characteristic  number  such 

that  for  any  other  characteristic  number  X. ,  <  lk| ,  [  5;  Chapter  6] . 
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A  sufficient  condition  for  X.  <  1  is  the  existence  of  a  nonnegative 
function  f(x)  such  that 

3.9  f(x)  <  /  K(x,  y)f(y)dy 

I 

for  x  in  I,  [  5;  Chapter  6]  . 

The  next  property  is. a  result  of  the  asymptotic  behavior  of  the 
n*  iterate  of  the  operator  K  when  X.Q  is  less  than  1. 

If  <  1 ,  then  for  any  nontrivial,  nonnegative  continuous  function  f(x)  there 
is  some  point  xQ  in  I  such  that 

3.10  f(xQ)  <  /k(xq,  y)f(y)dy  . 

We  eliminate  those  processes  having  no  absorption  by  assuming  a  positive 
probability  for  a  particle  initially  at  x  to  be  eventually  absorbed  at  T  or  in  IQ  , 

00 

3.11  0  <  a(x)  +  /{  /  b.(y,  s)p(x,  y;  s)g(s)ds  }dy 

I  0 


for  x  in  I . 
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We  can  now  state  the 

Theorem  3.  1.  Suppose  the  condition  3.  11  is  satisfied  in  addition  to.  the  conditions 
listed  in  Theorem  2.  1.  A  necessary  and  sufficient  condition  for  3.  6  to  have  a 
positive  continuous  solution  g(x)  different  from  the  function  l(x)  =  1  is 

3.12  XQ<1, 

where  is  the  unique  characteristic  number  determined  by  3.  8. 

We  apply  the  result  of  Theorem  3.  1  to  the  probability  fQ(x)  for  eventual 
extinction  defined  in  3.5  and  obtain  the 

Theorem  3.2.  Assuming  the  conditions  in  Theorem  3.1, 

(a)  if  X.g  >  1 ,  fQ(x)  =  1  for  x  in  I  and 

(b)  if  X  <  1  ,  fQ(x)  is  the  minimum  positive  solution  g(x)  to  3.6. 

Proof  for  Theorem  3.  1.  We  list  some  properties  for  the  operator  U  which  are  direct 
results  of  the  strict  monotonicity  of  h(y,  s;  t)  as  a  function  of  t  on  0  ' <  t  <  1 

and  the  normalization  1.7.  Letting  J  =  ,{x  in  -  il  a(x)  =  1  }  and  I  -  J=  {x  in  l|a(x)  <  1 }, 

we  then' have  the  following, for  any  nonnegative  continuous  function  f(x)’  which  is 
bounded  above  by  1  and  not  identically  1  , 
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'3.13  Uf(x)  =  1  for  x  in  J  and 

Uf ( x)  <  1  for  x  in  I  -  J . 

Suppose  g(x)  is  a  solution  to  3.  6,  g(x)  =  Ug(x)  .  Then  3.  13  implies 
1  -  g(x)  is  positive  for  x  in  I  -  J.  Since  l(x)  =  Ul(x)  ,  1.7,  1.11  and  2.  1 
imply 

00 

3.  14  0  <  1  -g(x)=  f  {  J  [  1  -h[y,  s;g(y)]]p(x,  y;  s)g(  s)  ds}  dy 

I  0 

</[l  -g(y)]K(x,  y)  dy 
I 

for  x  in  I  -  J.  Therefore  by  3.9,  the  characteristic  number  X  for  K(x,  y), 
determined  by  3.  8,  satisfies  <  1 ,  proving  the  necessity  of  the  condition. 

To  prove  the  sufficiency,  we  define  g^(x)  -  0  and  ' 

3.  15  g  (x)  =  Ug  (x) 

n+ 1  n 

for  x  in  I.  -  We  have  9.n(x)  -  1  for  x  in  J,  3,13.  Using  the  strict  montonicity 
of  h(y,  s;  T)  in  r  on  0  <  t  <  1,  we  have 
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3.16  Uf(x)  <  Ug(x) 

for  x  in  .  I  -  J  when  f(x)  <g(x)  on  I  -  J.  Assumption  3.  11  and  property  3.  13 
imply  0  <  g  ( x)  <1  on  I  -  J  ;  and  so,  3.  16  implies  g^x)  =  UgQ(x)  <  Ug^x)  =  g^(x) 
on  I  -  J .  By  induction  on  n ,  we  have  for  n  =  1 ,  2 ,  3 , 

3.  17  0  <g  (x)  <  g  (x)  <  1 

n  n+ 1 

for  x  in  I  -  J . 

The  inequality 


I  g  (x.)  -g  (x  )  1  <  I  a(x  }  -  a(x  )  1+  /Ik(x  ,  y)  -  K(x  ,  y)  I  dy 

n  1  n  2  —  1  1  t. 

and  the  continuity  of  K(x,  y)  show  g  ( x)  is  continuous  in  x  uniformly  in  n. 

n 

This  and  3.  17  imply  the  existence  of 

3.  18  lim  g  (x)  =  g(x) 

cl  -*■  oo  n 

uniformly  for  x  in  I.  Therefore  letting  n-*  00  in  3.  15  shows  g(x)  is  a 
continuous  solution  to  3.6  satisfying  0  <g(x)  <  1  for  x  in  I  -  J  and  g{x)  =  1 
for  x  in  J . 

We  now  show  g(x)  is  different  from  l(x)  =1  when  \  <  1  .  Tor  this 
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purpose  we  introduce  the  collection  of  kernels  defined  by 

00 

r  0 

K(x,y;r)  =  j  —  h(y,  s;t)p(x,  y;  s)g(  s)  ds 
0 

for  x,  y  in  I  and  0  <t  <  1  ;  so  that  K(x,  y;  1)  =  K(x,  y) ,  definedinl.il. 

We  use  the  assumptions  in  1.  1,  1.3,  1.4  and  1.7  to  state  K(x,  y;T)  is 
nonnegative  and  continuous  for  x,  y  in  I  and.  0<t<1.  Since  [j.(y,  s)  is 
positive  on  its  domain,  K(x,  y;T)  and  the  kernel 

oo 

P(x,  y)  =  /  p(x,  y;  s)  g(s)  ds 
0 

vanish  together. 

Each  kernel  K(x,  y;  t)  has  a  positive  characteristic  number  X  (t) 

O' 

determined  by  3.8.  Using  the  continuity  properties  of  K(x,  y;  t)  in  x,  y  and  t 
and  the  Fredholm  theory  for  integral  operators  with  continuous  kernels,  we  can 
show  that  the  Fredholm  determinant  d(X;T)  for  K(x,  y;T),  where  X  is  the 
spectral  parameter,  is  a  continuous  function  for  |  X.  |  <  °o  and  0  <  t  <  1  and  an 
entire  function  of  X  for  fixed  t.  This  is  sufficient  to  assert  Xq(t)  is  a 
continuous  function  on  0  <  t  <_  1 .  This  and  the  assumption  XQ  =  XQ(  1)  <  1 
imply  the  existence  of  a  T1  <  1  such  that  ^0(T)  <1  for  t  <  t  <  1.  We  now 
develop  a  contradiction  from  the  assumption  g(x)  =  l(x)  .  If  this  be  valid,  the 
uniform  convergence  of  Vjx)  to  g(x)  ,  3.18,  implies  the  existence  of  a  positive 
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integer  N(t  )  such  that  g  (x)  is  uniformly  close  to  l(x),g  (x)>  t  + - 1  , 

1  n  n  —  1  2 

for  n  >  N(t  )  .  Therefore  using  the  mean  value  theorem  and  the  strict  monotonicity  of 
a 

—  h(  y,  s;  t)  as  a  function  of  t  ,  we  have 

3. 19  (1  -  gN+1(x))  >  /(I  -  gN(y))  K(x,y;ri)  dy 

for  x  in  I  -  J.  Assuming  Xq(t^)  <  property  3.  10  gives  the  existence  of 
some  xQ  in  I  such  that 

/(l  -  gN(y))K(x0,y;T1)dy  >  1  -  g^C XQ)  • 

Since  K(x,  y;T^)  and  P( x,  y)  vanish  together,  we  know  xQ  is  in  I  -  J.  Therefore 
this  and  3.  19  imply 


WV 


>  i  - 


w 


or 


gN+l(V 


<3n(xo»' 


contradicting  3.17.  Consequently,  g(x)  is  different  from  l(x)  when  \  <  1, 
proving  the  sufficiency  of  the  condition. 

The  operators  U  and  T  have  a  structural  relation  which  is  utilized 
for  the  proof  of  Theorem  3.  2. 
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Proof  for  Theorem  3.  2.  If  X  >  1 ,  the  results  of  Theorem  3.  1  assert  f^(x)=  1, 
x  in  I,  proving  (a).  To  prove  (b)3  we  recall  the  probability  for  eventual  extinction 
is  given  by  f ^ ( >r)  =  lim  fQ(x,  t)  as  t  -*■  00  ,  3.5.  We  have  also  shown 
f(x,  t)  =limTgn(x3t)  as  n-*  00  with  the  operator  T  and  the  functions  g^(x,  t) 
defined  in  3.  3.  The  operator  T  has  a  monotone  property  similar  to  that  for  the 
operator  U,  3.1 6,  It  is 

3.  20  Tf(x ,  t)  <  Tg(x ,  t) 

for  x  in  I-J  and.  0  <t  if  f(x,t)  <g(x,t)  on  the  same  set.  Since  g  (x,  t)=0 
and  a(x3 1)  <  a(x)  3  1.  2,  we  have 

g^x,  t)  =  TgQ(x,  t)  <UgQ(x)  =  g^x) 

for  x  in  I-J  and  0  <  t  <  00 ,  and  so,  an  induction  argument  using  3.20  gives 

W1' l)  =  T9n,x’t)  tUg„(x>  =9„+i(x) 

on  the  same  set,  where  the  g  (x,  t)  are  given  by  3.  3  and  the  g  (x)  by  3.  15. 

n  n 


Therefore  we  have 
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3.21  fn(x,  t)  =  lim  g  (x,  t)  <  lim  g  (x)  =  g(x)  and 

u  ft  — >>co  n  —  ft. ->00  n 

f  (x)  =  lim  f  (x,  t)  <  g(x) 

0  00  o  — 

for  x  in  I  where  g(x)  is  defined  in  3.18. 

We  next  show  that  g(x)  is  the  minimal  positive  solution  to  3.6  . 

Suppose  h(x)  is  a  continuous  solution  to  3.6  satisfying  0  <  h(x)  <  1  for 
x  in  I.  Using  the  monotone  property  3.16  for  U,  we  have 

gx(x)  =  UgQ(x)  <  Uh( x)  =  h(x) 

for  x  in  I;  and  sousing  induction,  for  each  n 

g  (x)  =  Ug  (x)  <  Uh(x)  =  h(x) 
n+i  n  — 

for  x  in  I.  This  result  shows  lim  =  g  (x)  =  g(x)  <  h(x)  .  Therefore  3.21  implies 

n-^°o  n  — 

fQ(x)  =  g(x)  on  I,  completing  the  characterization  of  fQ(x)  as  the  minimal 
positive  solution  to  3.  6. 
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